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THE WORK REPORTED ON HEREIN was performed In 
support of the U.S. Department of Energy (DOE) 
Stirling Engine Highway Vehicle System Program 
for which the Li'wls Research Center (LeRC), 
through an Interagency agreement with DOE, has 
project management responsibility under the 
programmatic direction of the DOE Division of 
Transportation Energy Conservation. 

One objective of this program Is to make 
validated Stirling engine performance and de> 
sign optimization computer codes available to 
all who wish to do research and development 
work on Stirling engines. Although private 
Industry has been given contractual responsi- 
bility for the code development, NASA Lewis 
has been developing a performance code and 
attempting to validate It against experimental 
data. This approach makes It possible for us 
to contribute directly toward the development 
of validated codes and to manage the contrac- 
tual code development more Intelligently. 

One version of the code which models the 
engine of a rhombic-drive ground power unit 
(GPU) has already been reported in reference 
(1).* This paper updates that model and com- 
pares Its predictions with recent GPU test 
resu 1 ts . 

This test GPU was designed and built for 
the Army by General Motors. Tl>c unit was de- 
signed to produce 3 kW of electric power. 

GENERAL DESCRIPTION OF MODEL 

The GPU engine, for which this model pre- 
dicts performance, Is shown schematically In 
Fig. 1. The engine incorporates eight regen- 
erator/cocler units. The model assumes that 
Identical flow conditions exist In each of 
these eight parallel flow paths. The model 
represents the working space by a series of 
subdivisions called control volumes; this type 
of model Is sometimes called a nodal model. 

The engine working space Is represented 
by 13 control volumes, as shown In Fig. 2; 
the adjacent metal walls are represented by 
13 corresponding control volumes. The metal 
temperatures, except for those in the regen- 
erator, are assumed to be constant. This Is 
a reasonable assumption for any given run 
since the heater and cooler metal temperatures 


*Numbers in parentheses designate References 
at end of paper. 


are essentially boundary temperatures that are 
controlled by the combustor and the cooling 
water flow, respectively; these temperatures 
vary little over a cycle because the metal 
heat capacity is much greater than that of the 
gas. The appropriate fixed cooler metal tem- 
peratures are determined irom cooling water 
inlet temperature, flow rate and the calcu- 
lated heat rejection per cycle (this requires 
several iterations since heat rejection is a 
function of the cooler tube temperature); com- 
pression space metal temperature was assumed 
equal to the cooler metal temperatures. 

The model calculates Indicated power and 
efficiency for a given engine speed, mean 
pressure, heater and expansion space metal 
temperatures, and cooler water inlet tempera- 
ture and flow rate. The indicated efficiency 
is based on heat into the gas plus conduction 
losses. The model also simulates temperature, 
pressure, and flow variations over the cycle 
at various stations in the working space. 

The present model incorporates engine 
mechanical losses, as calculated from heat 
balances measured during the recent LeKC GPU 
tests . 

The thermodynamic equations and assump- 
tions used in developing the model were re- 
ported in reference 1. The modifications 
which have been made to the model since pub- 
lication of reference 1 are discussed in the 
appendix to this paper. 

The calculation procedure used in the 
model is outlined in Fig. 3. Each set of 
calculations, indicated in Fig. 3 within the 
inner loop (except for pressure-drop, conduc- 
tion, and shuttle losses), is made at each 
integration time step during each cycle. It 
is necessary to make the pressure drop calcu- 
lations only over the last cycle since, as 
explained in reference 1, the pressure drop 
calculations are decoupled from the heat and 
mass transfer calculations. Conduction and 
shuttle losses are calculated Just once during 
the last cycle. ( (Shuttle heat transfer oc- 
curs by heating of the displacer at the hot 
end of the stroke and cooling at the cold end 
of the stroke (2).) The regenerator metal 
temperature convergence scheme (1) was turned 
on after the Sth cycle, and off after the 25th 
cycle. During this time, the regenerator 
metal temperatures were corrected after each 
cycle and the cooler tube temperatures were 
corrected after every 4th cycle. With the 
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temperatures converged, die calculations pro- 
ceed through five mon cycles to allow further 
settling out and the run Is terminated with 
the completion of the 30th cycle. 

Thi.' computer program Is written In FOR- 
TR/\N V and, In card format. Is about 1300 
cards long (Including plotting subroutines). 
Current computing time Is about 3.73 minutes 
for 30 cycles on a L'NIVAC 1100, or 0.123 min- 
ute per cycle. This is based on 1000 Itera- 
tions per cycle (which Implies a time Increment 
of 2vl0"^ sec wlien the engine frequency Is 
30 Hz). 

The computing time can be decreased by 
decreasing the number of Iterations per cycle; 
the corresponding effect on predicted Indi- 
cated power and efficiency Is shown In Figs. 

4 and 3 for one senes of runs of an earlier 
version of the model. These figures show that 
some error was Introduced by reducing from 
1000 to 300 iterations per cycle and that the 
error became more significant if only 200 
Iterations per cycle were used. 

CM>ARISO.N' WITH OTHER MODELS 

The model described herein Is less gen- 
eral than the models of I'rlell (3,4), Shock 
(3), and Flnegold and Vanderbrug (6) In that 
It neglects the effects of gas Inertia (and 
pressure-wave dynamics) and gas kinetic energy. 
By neglecting these effects It is possible to 
use an Integration technique wlilch allows a 
trade-off between computing time and accuracy 
of solution without concern for numerical in- 
stability and which results In a savings of 
computer time. 

I'rlell has used his model to predict the 
performance of one prototype engine for which 
the above mentioned effects appear to be very 
significant at 30 Hz operating speed (7). 
However, these effects were assumed to be neg- 
ligible in the GPU (maximum Mach numbers were 
calculated to be about 0.03 for hydrogen and 
0.075 for helium at design spedd and pressure 
level (50 Hz, 6.90 10^ N/m2 or 3000 rpm, 

1000 psi)). 

The LeRC model differs from the Rios (8) 
and the Martini (2,9) models in that in more 
closely represents the "distributed parameter" 
nature of the working space by dividing each 
heat exchanger into several control volumes 
and making the heat exchanger inefficiencies 
an integral part of the cycle calculations. 
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The Lewis model is thus more general but less 
efficient, in terms of computing time, than 
the models of Rios and Martini. 

Nci tiler the LeRC or any one of the more 
general models (i.e., those of Urieii, Shock, 
and Fincgold and Vanderbrug) would be suitable 
for a design optimization program because of 
excessive computing time. One run of a 
Stirling engine design optimization program 
may require the prediction of engine perfor- 
mance at several thousand different configura- 
tions and/or operating conditions. Of the 
models mentioned here only the Martini or Rios 
models appear to be sufficiently efficient in 
terms of computing time to serve as the basis 
for an optimization code. 

ENGINE DESCRIPTION AND BACKGRaND 

The GPU-3 (Ground Power Unit) was ob- 
tained from the U.S. Army Mobility Equipment 
Research and Development Center (MERDC) at 
Fort Belvoir, Virginia. It was designed and 
built by the General Motors Research Labora- 
tories in 1965 for the U.S. Army. The prog- 
ress reports on the development of the engine 
and its components by General Motors as well 
as other useful information concerning the 
GPU are given in reference 10. 

The GPU-3 engine is a single-cylinder 
displacer- type engine with rhombic drive and 
sliding rod seals. It is capable of producing 
a maximum engine output of approximately 
7.5 kW (10 hp) with hydrogen working fluid at 
6.9x10^ N/m*^ (1000 psl) mean compression-space 
pressure. The piston swept volume is 120 cc 
(7.3 ln3). 

The engine obtained from Fort Belvoir 
was initially torn down and restored to oper- 
ating condition. It was then tested as part 
of the original GPU-3 package with only those 
changes made that were necessary to make the 
unit operable. Those results and a descrip- 
tion of the engine hardware are given in ref- 
erence 11. 

To obtain the data necessary for compari- 
son to the computer simulation, the engine was 
modified as follows: 

(1) The engine-driven accessories from 
the original GPU-3 package were removed with 
the exception of the lubrication oil system. 

(2) Air, water, fuel, and working fluid 
were provided from the facility support sys- 
tems . 
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(3) InsCrumontatlun was added to obtain 
an energy balance, engine temperature profiles 
working space gas temperatures, dynamic pres* 
sures, and to measure indicated work. 

The original CPU- 3 alternator and a sepa- 
rate resistance load bank were used to absorb 
the engine output power. The alternator was 
calibrated to define its efficiency at various 
speeds and output voltages. Since the GPU-3 
package was designed for a i-kW output, the 
alternator was not capable of absorbing the 
maximum output power of the engine. 

Engine dimensions and parameters needed 
in the simulation were determined and arc 
given in Table 1. A breakdown of the engine 
dead volume as calculated from hardware mea- 
surements is given in Table 2 together with 
the measured value of minimum working space 
volume (obtained by a gas displacement meth- 
od). Hie calculated value of minimum working 
space volume is seen to be within 2 percent 
of Che measured value. 

Tlie GPU-englne is shown schematically in 
Fig. 1. A more accurate schematic of the top 
portion of the engine (which shows the dimen- 
sions needed for heat conduction calculations) 
is shown in Fig. 6. The geometric relations 
between piston positions and crankshaft angle 
are shown schematically in Fig. 7. Equations 
for calculating compression-, expansion-, and 
buffer-space volumes are given in reference 1. 

PREDICTION PROCEDl’RE - The following pro- 
cedure was used in setting Che model up Co 
make comparisons with the experimental data: 

Heater Cube and expansion space metal 
temperatures were set to average values of 
the corresponding measured metal temperatures. 
Measured temperature profiles in the cylinder, 
regenerator, and insulation container (which 
encloses the insulation packed around the re- 
generators) were input to the model for use 
in conduction calculations. Measured water 
temperature into the cooler and water flow 
rate were used together with the calculated 
heat out to estimate the average cooler tube 
temperature (this required several iterations) 

The reference heat transfer and pressure 
drop CO relations identified in the appendix 
were used. However, the regenerator friction 
factor was adjusted to give best agreement 
between predicted and experimental Indicated 
powers. (This will be discussed further in 
the next section.) 
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Hie experlmencel ly determined mechanical 
losses, discussed in the appendix, were used 
to determine predicted brake power from pre- 
dicted indicated power. 

With the above procedure it was found 
that the time averaxed predicted gas tempera- 
tures were lower in the heater and higher in 
the cooler than the corresponding measured 
gas temperatures. Good agreement between pre- 
dicted and expi rimenta 1 average gas t.;mpera- 
tures could be obtained by increasing the 
model heat transfer coefficients in the heater 
and cooler. But increasing these coefficients 
made the agreement between predicted and ex- 
perimental engine power and efficiency worse. 
Tlicreforc, for these runs, we did not attempt 
to bring the predicted and experimental gas 
temperatures into close agreement. 

An optional procedure, which has been 
considered but not yet tried, would be to as- 
sume that the measured heater gas temperature 
is correct and to set the model heater metal 
temperatures to the values required to match 
the measured gas temperature. At this time 
It Is not clear which of the two sets of tem- 
perature measurements (metal or gas) is the 
more correct. 

COMPARISON OF MODEL PERFORMANCE PREDICTION 
WITH EXPERIMENTAL DATA 

Baseline tests were run to map the en- 
gine perfonr.ance over a range of heater gas 
temperatures from 595° to 705° C (1100° to 
1300° F), mean compression-space pressures of 
1.4 to 6.9x10^ N/m^ (200 to 1000 psl), and 
engine speeds from 1000 to 3500 rpm with both 
hydrogen and helium as the working fluids. 
These experimental results are reported in 
reference 12. 

Comparison of experimental and predicted 
indicated powers suggested that the engine 
pressure drop was underestimated when the ref- 
erence set of friction factors and end effects 
were used in the model. (The end effects 
calculations are described in ref. 1.) Since 
the regenerator accounts for most of the cal- 
culated frictional pressure drop, the regen- 
erator friction factor was adjusted to get 
better agre -ment between predictions and data. 
For runs made with hydrogen working gas the 
plot of indicated power shown in Fig. 8 sug- 
gests that predicted indicated power agrees 
best with the test data if four times the 
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'’sference regenerator friction factor it uaed; 
thla factor of 4 produced the best overall 
agreement for all tc-st points made with hydro* 
gen working gas. 

Tlie reference regenerator and cooler 
friction factors are based on steady-flow air 
tests that were made on Che newly manufactured 
regenerator-cooler units prior Co Che engine 
tests. Examination of Che regenerator-cooler 
units after the engine tests revealed they had 
been contaminated with oil. Steady flow air 
tests of the oil contaminated components 
(i.e., tests on (1) the complete heater head- 
regeneraCor-cooler flow system and (2) three 
individual regenerator-cooler units) were com- 
pleted after Che computer predictions dis- 
cussed in this paper were completed. Prelimi- 
nary analysis of these results indicate Chat 
the three contaminated regencraCcr-cooler com- 
ponents have about one-third greater pressure 
drop than wtien clean and also the spread in 
pressure drop from one unit to another has 
increased. Also the flow tests on Che com- 
plete heater head, regenerator, cooler flow 
system suggest that the pressure drop due Co 
end effects appears almost as large as the 
regenerator frictional pressure drop. It 
thus appears that loss due Co end effects is 
being underestimated in the computer model. 
Further analysis using these airflow results 
will be required before the appropriate modi- 
fications can be made to the computer model. 

ocher possible explanations for addi- 
tional pressure drop are; 

(1) Increased pressure drop due to 
methane (which can be produced by Chennai 
cracking of the oil contaminant). In refer- 
ence 13 (p. 9) it is stated Chat in the Gen- 
eral Motor tests "the effect [of methane 
production] was much more significant than 
had ever been noted by plugging of the heat 
exchangers by solid residue." In future 
tests it is planned to periodically analyze 
samples of the working gas to evaluate the 
effect of such impurities. 

(2) Error due to use of steady flow 
friction factor correlations to predict pres- 
sure drop for periodically reversing flow. 

(3) Errors in pressure drop prediction 
due to the neglect of gas inertia and gas 
kinetic energy. (These effects are probably 
not significant in the GPU at the low Mach 
numbers encountered (< 0.1) when hydrogen 


7 



and helium are used at apeeda at or below 
3500 rpm.) 

Another poaalbllity la that we may be 
undurestimatlnK the true mechanical loai; 
auch an error would make the apparent prcaaurc 
drop loaa appear larger (ace the diacuaaion of 
mechanical loaa in the appendix). 

Performance predicCiona are compared with 
the experimental reaulta for hydrogen working 
gas in Figa. 9 to 12; four times the reference 
regenerator friction factor waa used in the 
model to make all predictions shown in these 
figures. Brake power is shown as a function 
of engine speed for several pressure levels 
in Fig. 9. The 4.14x10^ N/m^ (faOO pal) pres- 
sure was the highest experimental pressure 
that could be achieved without overloading the 
alternator; the alternator limitation also 
prevented completion of the experimental curve 
at this highest pressure level (sec Fig. 9). 
Agreement is within 6 percent of the experi- 
mental values at the 2.7bvlOb N/m2 (400 psi) 
and 4.14x10*’ N/m^ (600 psi) pressures. At 
the 1.38xl0*> N/m2 (200 psi) mean pressure 
level although the absolute error In predicted 
power is no greater than 0.22 kW (0.3 hp) the 
percentage error is as high as 44 percent due 
to the low power level. 

Predicted and experimental brake effi- 
ciencies arc compared in Fig. 10 for the same 
hydrogen runs. Brake efficiency is defined 
as follows: 

Brake efficiency 

brake power 

rate of heat input to engine 

Two sets of experimental brake efficiencies 
are shown for each pressure level. Results 
of two methods of determining the experimental 
heat input rate, used in determining these two 
sets of experimental brake efficiencies, are 
compared with the predicted heat input rate 
in Fig. 11. One set of values was determined 
by a hot end energy balance as follows: 

Heat input ratel. ^ , 

'hot end energy balance 

> heat produced by fuel 

- [sum of heat loss rates via ex- 
haust gas, fuel nozzle cooling 
water, conduction, and radiation] 
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Till* leC of heat Input valuea waa the higher 
and therefore gave the l'>wer experimental ef* 
ficleticiea. Complete combustion of the fuel 
was assumed. Hie lower hejt Input values In 
Fig. 11, wlilch correspond to the higher ex* 
pcrimental efficiencies in Fig. 12, were de- 
termined as follows: 

h‘*.it Input rate I ... . , 

'cold end energy balance 

■ brake power [sum of heat rejection 
via the cooling water and oil] 

In Fig. 10 the predicted value of effi- 
ciency at ^.lAxlO*’ N'/m^ (600 psl) Is 7 percent 
above the upper experimental value. At 
2.76x10^ N/m^ (AOO psl) the predicted values 
of efficiency range from 6 to lb percent 
above the upper experimental values with best 
agreement at the higher speeds where power 
and heat input are largest. At 1.38x10^ N/m^ 
(200 psl) there Is a large spread between the 
two sets of experimental va.ues Indicating a 
high degree of uncertainty. This Is probably 
due to Che fact that Che errors In estimating 
various losses (radiation, convection, exhaust 
losses, etc.) represent a larger percentage 
of Che brake power and heat Input to Che en- 
gine than at higher pressure levels. The pre- 
dicted vnlues of efficiency range from 6 to 
30 percent above the upper set of experimental 
values at this mean pressure with best agree- 
ment again at Che higher speeds. 

In Fig. 11 the difference between ex- 
perimental and predicted values of heat input 
may be attributable In part to lowered exper- 
imental regenerator effectiveness (and there- 
fore higher heat input) due Co oil contamina- 
tion. At 4.14x10^ N/m^ (600 psl) the pre- 
dicted heat Input is about 3 percent below 
Che lower experimental value. AC 2.76x10^ 

N/m^ (400 psl) the predicted hecc Input Is 
about 0.43 kU (0.6 hp) lower over Che whole 
range. At 1.38x10^ N/m2 (200 psi) the pre- 
dicted heat Input is within Che range of the 
experimental values except in the low speed 
range where it falls slightly below the lower 
experimental values. 

The effect of different heater gas tem- 
peratures on Che comparison of predicted and 
experimental brake powers is shown in Fig. 12 
for hydrogen working gas at 2.76x10^ N/m^ 

(400 psi). The comparison suggests Chat the 
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model accurately predicted the effects of 
chanKca In heater gaa temperature on brake 
power. Correlation of predicted and experi* 
mental values of brake efficiency and heat ln> 
put tor h49® C (1200” F) and 593" C (1100° F) 
were nut significantly different than for the 
704° C (1300° F) heater gas temperature. 

For helium working gas, comparisons of 
predicted and experimental Indicated power 
such as that shown In Fig. 13 show that best 
agreement with the test data results when 
about 2.6 times the reference generator fric- 
tion factor Is used. Although the comparison 
Is shown for 2.76>10" N/m^ (400 psl) mean 
pressure and 704° C (1300° F) heater gas tem- 
perature, the choice of the factor of 2.6 was 
based on such comparisons over the mean pres- 
sure range from 2. 76^10” N/m^ (400 psl) to 
6.9(M0*» N/m2 (1000 psl) at 649° C (1200° F) 
heater gas temperature. 

Theoretically the same friction factor - 
Reynolds number correlations and end effects 
calculations should apply to both hydrogen 
and helium. The fact that different "correc- 
tion factors" for regenerator friction factor 
were required to get best agreement with the 
data may be due In part to the following: 

(1) rhe hydrogen tests were made after 
the helium tests; thus the larger friction 
factor required to fit the hydrogen data may 
be due In part to additional solid residue 
from oil that was deposited after completion 
of the helium tests. 

(2) Also any methane produced by break- 
down of the oil would cause a larger percent- 
age error In the pressure drop calculation 
for the hydrogen than for the helium runs 
(since methane Is closer In density to helium 
than to hydrogen). 

(3) It Is shown In the appendix that the 
experimentally determined mechanical losses 
were higher for helium than for hydrogen. 

This result Is supported by results discussed 
briefly In reference 10 (Sec. 1.206, p. 4). 

If, however, the results shown In the appendix 
overestimate the difference between helium 
and hydrogen mechanical losses, an effect 
would be to cause an apparent difference be- 
tween friction factors for hydrogen and helium 
such as chat already noted. 

Performance comparisons for helium work- 
ing gas are shown In Figs. l4 to 18; 2.6 times 
the reference regenerator friction factor was 
used in the model to make all predictions 
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•ho%m In thcav figure*. In Che br*Ke power 
plot of Fig. 14 (for 649° C (1200° F) he*ter 
temperature) the predicted power range* from 
2 to 15 percent above the exper imental power 
for all point* except the one of lowe*t power. 

The correaponding brake cfficiencie* are 
compared in Fig. 15. In Fig. 15(a) the pre* 
dieted effiriencie* range from 3 to 17 percent 
above the C|:p.' r experimental value* at all 
points except the point of lowest efficiency. 

In Fig. 15(b) for 5.52vlO<> N/m^ (800 pal) 
the predicted value* lie above the tipper ex* 
perimental value* by from 10 to 15 percent 
over the range tested. In Fig. 15(c) for 
6. 90^:10** N/m* (1000 pal) the predicted values 
exceed the upper experimental value* by about 
12 percent over the range tested. 

The corresponding value* of heat in are 
compared in Fig. 16. At 6.9(N10*’ N/m2 (1000 
p*l) and 5.52vlO^ N/m^ (BOO p*l) the predicted 
values lie within the ranee of experimental 
value*. For 4.l4vl0^ N/m^ (600 pal) and 
2.76x10*^ N/m2 (400 psl) the predicted heat in 
is slightly low at the higher engine speed* 
but still within 4 percent of the lower ,'X> 
perimental values. Thus errors in predicted 
brake efficiency for this set of helium runs 
are primarily due to overprediction of power. 

Predicted and experimental brake powers 
.‘or 704° C (1300° F) nominal heater gas tem- 
perature are compared io Fig. 17. The shapes 
of the predicted curves are in good agreement 
with the data; however, the "''onitudes are 
more optimistic than for those runs made at 
649° C (1200° F) nominal heater gas tempera- 
ture. At 5.52x10^ N/m2 (800 psl) the pre- 
dicted values are within 10 percent of the 
experimental values. However at 4.14x10^ 

N/m^ (600 psl) and 2.76x10° N/m^ (400 psl) 
the predicted values range from 19 to 28 per- 
cent and from 7 to 34 percent, respectively, 
above the experimental values. Correlation 
between predicted and experimental power at 
4.14x10^ N/m2 (600 psi) is notably worse than 
for the corresponding 649° C (1200° F) case 
shown In Fig. 14. At 1.38x10^ N/m2 (200 psl), 
the maximum absolute error is 0.15 kW (0.2 hp) 
over the range tested. 

Predicted and experimental brake powers 
for 593° C (1100° F) nominal heater gas tem- 
perature are compared in Fig. 18. 


11 





Sl'MMARY 


Tills paper compjres predictions of an up- 
dated version of the I^RC GPl' Stirling engine 
computer model (1) with recent GPU tests per- 
formed at LcRC. 

Comparisons ol predicted and experimental 
Indicated powers Implied that the reference 
pressure drop calculations were underestimat- 
ing actual engine pressure drop. This result 
seemed consistent with (1) the presence of oil 
contamination In the regenerator-cooler units 
after engine testing, (2) preliminary analysis 
of flow tests on the entire flow system and on 
several of the oil contaminated regenerator- 
cooler units, and (3) the possible presence 
of methane In the working pas due to thcm.al 
cracking of oil. 

Another possibility Is that the LeRC 
mechanical loss data may underestimate the 
true mechanical loss; such an ?rror would 
make the apparent pressure drop loss, as im- 
plied by exp'^ 1 . imental Indicated power, appear 
larger . 

The procedure used In this comparison 
was to adjust the regenerator friction factor 
to give best agreement between predicted and 
Indicated powers. For hydrogen working gas 
a factor ot 4 times the reference regenerator 
friction gave best results and for helium a 
factor of about 2.6 was best. Tills apparent 
difference In friction factors for hydrogen 
ana helium, although theoretically unreason- 
able, would be consistent with (1) progressive 
accumulation of solid residue from oil con- 
tamination (the hydrogen tests were done after 
the helium tests) and (2) possible presence 
of methane (wtilch would affect the hydiogen 
results more than the helium. Any overesti- 
mate of the true difference between helium 
and hydrogen mechanical losses would also 
tend to produce an apparent difference In he- 
lium and hydrogen friction factors such as 
that noted above. 

Results of the comparisons between pre- 
dicted and experimental brake power and effi- 
ciency are sumiarlzcd In Table 3; the predic- 
tions were made using the adjustments In re- 
generator friction factor discussed above. 
Table 3 shows that for hydrogen working gas 
the predicted values of brake power ranged 
from 0 to 6 percent above the experimental 
values at the two higher pressure levels. 

For the corresponding predicted orake effl- 
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ciencies the error range was higher (from 
6 Co 16 percent above the upper set of expert* 
mental brake efficiencies). In general the 
highest percentage errors occurred at Che low 
power and efficiency points. For helium with 
649° C (1200° F) heater gas temperature the 
predicted brake power was 2 to 15 percent 
greeter than the experimental (except for one 
point) over a wider pressure range Chan could 
be covered with hydrogen working gas; Che per* 
cent error range for the corresponding pre* 
dieted brake efficiencies was about Che same 
as for brake power. The correlation between 
predicted and experimental brake powers is 
seen, in Table 3, to be worse for helium at 
704° C (1300° F) heater gas temperature than 
at 649° C (1200° F). 

Due to uncertainties in Che various true 
losses * pressure drop, mechanical loss, leak* 
age, heat exchanger inefficiencies, etc. * it 
is possible to obtain good performance predic* 
Cions for a given set of engine tests with a 
combination of incorrect loss predictions. 

Thus individual losses are better determined 
by special tests and/or instrumentation rather 
Chan implied by overall results. TIterefore 
attempts must be made to obtain better defini* 
Cion of Che various losses. Tlie recent flow 
tests on the oil contaminated components imply 
that pressure drop loss due to end effects and 
oil contamination has not been properly ac- 
counted for in the model. Periodic samples 
of Che working gas must be made to determine 
if significant amounts of methane or other 
contaminants are present. Motoring tests of 
the engine will be made and an attempt will 
be made to Improve the direct measurement of 
indicated power to develop more confidence in 
the mechanical loss definition. It may be 
possible to improve the leakage lose estimate 
by use of techniques discussed in reference 10 
(Sec. 5). 

It should be noted that only one of the 
test points discussed in this paper was taken 
at the design operating pressure of 6.90x10^ 
N/m^ (1000 psi). It may be unreasonable to 
expect great accuracy at pressures as low as 
1.38x10° N/m^ (200 psi) when several tenths 
of a kilowatt can lead to such large percent- 
age errors. The next series of tests with 
Che engine on a dynamometer should result in 
complete maps of engine performance for use 
in model validation. 
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APPLNDIX - MODEL UPDATES SINCE 
PUBLICATION OF REFERENCE I 


1. CALCULATION OF COOLER H'BE TEMPERA- 
TURE - When the Model was reported In refer- 
ence I It w<<i> assumed that the water side 
heat transfer coefficient in the cooler was 
so large compared to the gas side coefficient 
that the cooler tube temperature was essenti- 
ally at Che average water temperature. Infor- 
mation given in reference 10 suggests Chat the 
above assumption was a poor one. Consequently, 
cooler Cube temperature is now calculated as 
follows : 

The heat transfer correlation used on 
the water side is 


k 


. ^0.55,.0.33 

*^Re Sr 


( 1 ) 


(from (15)) wliere 


h water side heat transfer coefficient 

w 

D hydraulic diameter (lieater Cube out- 

side diameter) 

k water thermal conductivity 

N' _ Reynolds number 

Kb 

N Prandtl number 

I K 

In addition a fouling factor F, of 0.882 
ra2-°C/kW (0.005 hr-f t^-^R/Btu) is assumed to 
exist on Che water side. Thus Che total co- 
efficient on Che water side is; 


‘'w 


Thus if h,^ - 1.134 kW/m2-°C (2000 Btu/ 
hr-ft^-°R) then 


1 

1.134 


J 

+ 0.882 


0.567 kW/m^-°C 


(1000 Btu/hr-ft2-°R) 

The fouling factor, in this case, cuts Che 
water side coefficient in half. In calcu- 
lating Che water side Reynolds number 
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the flow per unit area Ci la based on an ef- 
fective flow area of 3.01 cin^ (0.467 In^) tor 
the total cooler water flow (one-half of the 
cooler water flow passes through each cooler). 

The equation used to calculate tube tem- 
perature from the Inlet water temperature is 

^tube * 


+ Q, 


out 


2CpwWw 


ln(rQ/r^) 

+ — . _ . + 


2i: Lk, 


m 



( 2 ) 


where 




Qout 

heat out of engine 

over one 

eye le 

^PW 

specific heat ot water 


Ww 

water flow rate 




outside and inside 
tube 

radii of 

cooler 


L total heat transfer length of cooler 

tubes 

kjj^ thermal conductivity of cooler tube 

h(. total heat transfer coefficient, 

water side 

total heat transfer area on water 
side 


Since the value of is changing 

from cycle to cycle during the convergence 
procedure it is necessary to repeat this cal- 
culation several times. The correction is 
made after every 4th cycle during the se- 
quence of cycles when the regenerator metal 
temperatures are being corrected after each 
complete cycle. 

The portion of the heat rejection due to 
pressure loss was neglected in determining 
the temperature difference between the cooler 
wall and the cooling water. This caused an 
overprediction of power estimated to be no 
more than 2 to 3 percent (with the largest 
errors occurring with helium working gas at 
the highest engine speed). 

2. LEAKAGE LOSS - Tlie simplified equa- 
tions incorporated in the model for calcu- 
lating leakage throughout the cycle, de- 
scribed in reference 1, were not used for 
this set of runs. These equations always 
give the same mean pressure for working and 
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buffer spaces. Ttte test data showed that this 
assumption Is not a good one. 

Instead an estimate of the power loss due 
to leakage between working and buffer spaces 
wij made based on Information given in refer- 
ence 10 (Sec. 6.006). Accordin.*, to this In- 
formation the power loss due to leakage is 
about 0.43 kW (O.P hp) for hydrogen working 
gas at the design mean pressure of 6.90^10^ 

N/m^ (1000 psi). An equation derived in this 
reference for power loss due to leakage be- 
tween working and buffer spaces implies that 
this power loss is proportional to the square 
of the mean pressure level and is independent 
of engine speed. The equation also indicates 
that this leakage power loss is inversely pro- 
portional to the square root oi the molecular 
weight. Using these relationships and assum- 
ing that the 0.43 kW (0.6 hp) leakage loss for 
hydrogen is due only to leakage between work- 
ing and buffer spaces, the leakage losses in 
Table 4 were generated and used in the model. 

3. MECHANICAL LOSS - Mechanical loss 
curves for hydrogen and helium are shown in 
Figs. 19 and 20. At the 4.14x10^ N/m^ (600 
psl), 38.33 Hz (3300 rpm) point the mechanical 
loss with helium working gas is about 19 per- 
cent greater than with hydrogen. At 2.76x10^ 
N/m^ (400 psi) mean pressure the mechanical 
loss with helium is higher than with hydrogen 
by from 9 percent at 16.67 Hz (1000 rpm) to 
26 percent at 38.33 Hz (3300 rpm). At 
1.38x10^ N/m2 (200 psi) the mechanical loss 
with helium ranges from 7 to 42 percent higher 
than with hydrogen over the speed range tested. 

These mechanical loss curves were deter- 
mined from experimental heat balances taken 
during the recent LeRC GPU tests. It was as- 
sumed that the sum of the heat to Che buffer 
space cooling water and Co Che oil equaled 
the mechanical loss. According to refer- 
ence 10 (Sec. 6.006) General Motors concluded 
that this heat balance approach to estimating 
mechanical loss was probably more accurate 
Chan use of motoring loss data. 

The mechanical loss for hydrogen at 
2.76x10^ N/m.2 (400 psl), from Fig. 19, can be 
compared to General Motors motoring loss 
data (10) (Sec. 6.006, p. 36); at 3000 rpm 
the motoring loss data, which Includes leak- 
age, implies a loss about 0.37 kW (0.3 hp) 
greater than the mechanical loss data of 
Fig. 19. Table 4 shows a leakage loss of 
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0.07 kW (0.10 hp) for hydrogen at 2.76x10^ 

N/ra^ (400 psi). 

The implication of the above is that the 
LeRC mechanical loss data indicates about 
0.3 kW (0.4 hp) less mechanical loss (at 3000 
rpm and 2.76x10^ N/m^ (400 psi)) than the 
General Motors motoring loss data. The dif- 
ference between the mechanical losses pre- 
dicted by the two sets of data decreases lin- 
early with speed until they agree at approxi- 
mately 1600 rpm. 

The relative difference between the he- 
lium and hydrogen mechanical losses (Figs. 19 
and 20) agrees well with information given in 
reference 10 (Sec. 1.206, p. 4) for 3000 rpm 
and 4.83x10^ N/m^ (700 psi). 

4. HEAT TRANSFER CORRELATIONS - Steady- 
flow heat transfer correlations from refer- 
ence 14 are used in the heater and cooler 
(Fig. 7-1, p. 123 of ref. 14) and the regen- 
erator (Fig. 7-8, p. 129 of ref. 14). The 
validity of these correlations is more widely 
accepted than the ones used in the model when 
it was reported on in reference 1. 

5. FRICTION FACTOR CORRELATIONS - Curves 
of friction factor, f, as a function of Reyn- 
olds number, Nj^g, were derived from NASA Lewis 
steady-flow air data for the regenerators and 
coolers; these curves are shown in Figs. 21 
and 22. The regenerator Reynolds number is 
based on hydraulic diameter (rather than wire 
diameter). The friction factor correlation 
used in the heater is a standard one from 
reference 15; the form of this correlation 
used in the model is : 

f = 16/Nj^^ Nj^^ < 1500 

f = 0.046/N°g^ 1500 

This correlation was not satisfactory for use 
in the cooler because the small cooler tube 
diameter implies a large relative roughness. 
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T«bl« 1 • CPV 3*2 engine Diaensioos mod fmtmmrntmf 


Cylinder bore et liner, m (in.) b.tf (2.731) 

Cylinder bore above liner,* cb (in.) 7.01 (2.76) 

Cooler 

Tube length, cm (in.) ••••, 4.61 (1.013) 

Heat transfer length, cm (In.) 3.35 (1.366) 

Tube inside dimeter, cm (In.) 0.108 (0.0423) 

Tube cMcsidc diameter, cm (in.) 0.136 (0.0623) 

Number of tubes per cylinder 

(or number of tubes per regenerator) ••..•.•.•••. 312 (36) 

Heater 

Nran tube length, cm (in.) 24.33 (6.63S) 

Heat transfer length, cm (in.). •.•.•«•....••• 13.34 (6.12) 

Cylinder tube, cm (in.) 11.64 (4.303) 

legenerator tube, cm (in.) 12.86 (3.073) 

Tube Inside diameter, cm (in.) 0.302 (0.116) 

Tube outside diameter, cm (in.) ••• 0.403 (0.16) 

Nuaber of tubes per cylinder 

(or ntsiber of tubes per regenerator) .40 (3) 

Cold end connecting ducts 

Uf«th, cm (in.) 1.36 (0.625) 

Duct inside diameter, cm (in.) 0.367 (0.233) 

Nuaber of ducts per cylinder 0 

Cooler end cap, as’ (in^) 0.276 (0.0170) 

tegenerators 

langth (inside), cm (in.) 2.26 (0.89) 

Diameter (inside), cm (in.) 2.26 ^.89) 

Number per cylinder 0 

Material Stainless steel wire cloth 

Nusber of wires, per cm (per in.) 79*79 (2(KK200) 

Wire diameter, cm (in.) 0.004 (0.0016) 

Number of layers 308 

Filler factor, percent •••• 30.) 

Angle of rotation between adjacent screens, deg ...•••»••.. 3 

Drive 

Connecting rod length, cm (in.) 4.60 (1.810) 

Crank - radius, cm (in.) 1.38 (0.34)) 

Eccentricity, cm (In.) 2.08 (0.820) 

Miscel laneoos 

Displacer rod dimseter, cm (in.) 0.632 (0.373) 

Piston rod diameter, cm (in.) 2.22 (0.873) 

Displacer diameter, cm (in.) 6.96 (2.740) 

Displacer vail thickness, cm (in.) 0.139 (0.0625) 

Displaf.cr stroke, cm (in.) 3.12 (1.23) 

Expansion space clearance, cm (in.) 0.163 (0.064) 

Compression space clearance, cm (in.) • 0.0)0 (0.012) 

Buff.r .pKc Mxiui volw. (!.}) S21 (31.78) 

Total working space minimum voliae, cm (in )...... 23). 3 (14.23) 

*Top of displacer seal is at top of liner at displacer TDC. 


E- ;r^r 


T.bl. 2 - GK-i Stlrlti^ bflor D..4 Volu 
[Voliaes are given in cm cm (cm tn.)| 
space clearamce volume 


Displscer clearamce (aroumd dlsplecer) 

Clcarence voltmm above displacer 
Volume from end of heater tubes into cylinder 
Total 

11. Heater dead voliaee 

Insulated portion ol heater tubes next to 
expansion space 
Heated portion of heater t u bes 
Insuleted portion of heater tubes next to 
regenerator 

Additianel volume in (our iseeter tubes weed (or 
instriaentet ion 
Volums in header 
Total 

111. Rexenerstor dead volusir 

Entrance valimm into regeneretors 
Volume within matrix and retaining disks 
Volume between regenerators mod coolers 
Voltae la snap ring grooves et cmd of coolers 
Totel 

IV. Cooler dead wolxmm 

Volume in cooler tubes 
V. Cumprssslon in spece clearamce volume 

Exit voluam from cooler 

Volume in cooler end caps 

Volume in cold cod commectiag ducts 

Power piston cleerente (eround p o w er piston) 

Clearamce volume between displecer amd powe r pistom 

Votise at connect toms to cooler end caps 

Volume In piston "notches** 

Volume eround rod In botton of displacer 
Total 

Total dead voit»r 
Minimal live voliane 

Calculated sdniaM total working space 
Volume 


3.34 (0.204) 
7.41 (0.432) 
1.74 (0.104) 

T775 (Sr^*2) 


9.68 (0.361) 

47.46 (2.866) 
LJ.26 (0.811) 

2.74 (0.167) 

7.»7 (0.4P8) 
80.8 (^.61)) 


7.36 (0.446) 

33.4 (1.238) 

2.36 (0.138) 


Measured value of niniaam total working space volume 
(by vdlime displacement) 

Change ia working space voliase due to minor engine 
mediticacion 


63 

3 

(S.**8) 

13 

33 

(0.801) 

3 

*2 

(0.2») 

2. 

77 

(0.3*8) 

3. 

36 

(0.217) 

7, 

28 

(0..i3) 

1 

14 

(O.O.'O) 

2 

33 

(0.1*2) 

0 

06 

(0.004) 

0 

11 

(0.00.) 

21 

18 

(1.7*3) 

163 

13 

(11.787) 

1*1 

18 

(;.3*j) 

232. 

3 

(U.178) 

233 

3 

0*.23) 

2 

3 

(0.1» 


2M.0 (14.40) 



Table 3 - Sunnary of Cotaparisons Between Predicted and Experiaental Brake Power and Efficiency 


Working gas 

Heater 

gas 

temperature, 
“C (®F) 

Performance 

paraawter 

Mean pressure, 
N/n2 (pal) 

Percent or absolute value of predicted 
performance paraiaeter relative to 
experinental value 

Hydrogen 

704 (1300) 

Power 

4.14x10* (600) 

0 to 6 percent above 




and 





2.76x10*' (400) 





1.38xl0*> (200) 

0.07 kW (0.1 hp) to 0.22 kW (0.3 hp) above 



Efficiency 

4.14X10** (600) 

7 percent above upper experianntal value 




2.76x10*' (400) 

6 to 16 percent above upper experimental 





value 




1.38x10* (200) 

6 to 30 percent above upper experlawntal 





value 

Heliun 

649 (1200) 

Power 

6.90x10* (1000) 

2 to 15 percent above (except for lowest 




to 

pewer point) 




2.76x100 (400) 




Efficiency 

6.90x10* (1000) 

3 to 17 percent above upper experimental 




to 

value (except for lowest efficiency 




2.76x10* (400) 

point) 


704 (1300) 

Power 

5.52x10* (800) 

7 to 10 percent above 




4.14x10* (600) 

7 to 34 percent above 




and 





2.76x10* (400) 



Table 4 - Leakage Losses for Hydrogen and Heliun 


Hean pressure, 
N/n^ (psi) 

Hydrogen loss, 
kW (hp) 

Heliun loss, 
kW (hp) 

6.90x10* 

(1000) 

0.45 

(0.60) 

0.31 

(0.42) 

5.52x10* 

(800) 

0.28 

(0.38) 

0.20 

(0.27) 

4 . 14x 10* 

(600) 

0.16 

(0.22) 

0.12 

(0.16) 

2.76x10* 

(400) 

0.07 

(0. 10) 

0.05 

(0.07) 

l.38x 10* 

(200) 

0.015 

(0.02) 

0.007 

(O.Ol) 
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Figure 1. * Schematic of a single-cylinder Stirling engine 
with rhombic drive. 
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Figure 2 . - Heat- and mass-transfer control volumes. 


Inputv Frequency of optrjdon 

Pressure level «f initiel volume) 
InitMl temqerjtures 
Coolinq water flow rale 


Make 

revision 


Revise cooler lube 
lemperalure? 


19. Revise 

reqenerator metal 
Icmperatures 
to speed up 
convergence 


I. Tine • Time « Allimei 

Z update pressure level, P • f (volumes, temperatures, masses) 

3i Update temperatures in eKh gas control volume for effect of change 
in pressure level 

4 Update mass distribution 
Update How rates 

6. Update temperatures in eKh gas control volume lor efiKt of How 
between control volumes 

7. Update heal -transfer coeflicients between metal and gas control 
volumes 

8. update temperatures in eKh gas control volume lor etiKt of heat 
transfer ibelWMn gas and metal) 

9. Update regenerator metal temperatures for eflKt of heat transfer 
(between gas and metal) 

la Calculate pressure drops inKessary over only one cycle) 

II. Calculate conduction and shuttle losses i|ust once during last cycle) 

12. Sum up heat transfers between gas and metal (or eKh component 

II Update leakage How between working and buffer spKes loptionaH 

14. Update erpans ion-, compression*, and buHer-spKe volumes 

15. Calculate work done in .'vpansion and compression spaces and sum 
up to get total indicated work 


- 15. Is cycle complete? 

I Yes 

16. Calculate indicated power and efficiency 

17. Have spKilied number of cycles been completed? 


21. Calculate brake power and efficiency 
Figure i - Outline of calculation procedure. 
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NUMBER OF ITERAnONS PER CYCLE 

Figure 4. - Indicated power and change in indicated power as 
function of number of iterations per cycle. Working fluid, 
helium; frequency, $0 hertZ; heater-tube temperature. 

978 K (176DP RH cooler-hibe tmperature. 325 K i585° Rh 
mean pressure. 4.43xl(rNfm‘ i643psil. 



NUMBER OF ITERATIONS PER CYCLE 

Figure 5. * Indicated efficiency and change in indicated efficiency 
as a function of number of iterations per cycle. 
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Figure 7. ~ Schematic showing geometric relations between piston 
positions and crankshaft angle. 
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Figure 8. - Effect of regenerator friction factor, 
*REG ' ^ predKted indicated power. 
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Figure 9. - Brake power as a (unction of engine 
speed at several mean pressures. 
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Figure 10. - Brake efficiency as a 
function of engine speed at 
several mean pressures. 
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Figure 12. - Brake Dower as a function of speed 
at three different heater gas temperatures. 
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Figure 11. - Rate of heat input as a function of 
engine speed. 
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Figure 13. - Effect of regeneretor friction factor. 
fpEG . on predicted indicated power. 
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Figure 14. - Brake power as a function of en- 
gine speed at several mean pressures. 
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Figure 16. 
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Figure 17. * Brake power as a furKtion of engine 
speed at several mean pressures. 
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Figure 18. - Brake power as a function of en* 
gine speed at 593P C tllOO^ Fi nominal 
healer gas temperature. 
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Figure 20. * Mechanical loss as a function of engine 
speed for helium working gas. (Determined from 
experimental heat balance. ) 
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Figure PI. • Regenerator friction factor, f as a function of 
Reynolds number. Nr,. 



Figure 22 . * Cooler friction factor, f. as a function of Reynolds 
number. Nr,, 



